Intriguing examples of inhomogeneous broadening 
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Abstract 

Three problems are considered in which inhomogeneous broadening can yield unusual conse- 
quences. One problem involves the energy levels of atoms moving within nanopores of nearly 
cylindrical cross section. A second involves atomic or molecular motion in a quasi-one dimensional 
interstitial channel within a bundle of carbon nanotubes. The third problem involves motion within 
a groove between two nanotubes at the surface of such a bundle. In each case, the density of states 
at low energy is qualitatively different from that occurring in the perfectly homogeneous case. 
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1. INTRODUCTION 



Problems involving inhomogeneous broadening are ubiquitous in physics and chemistry. 
The usual consequence of particles experiencing a variety of environments (or a range of 
velocities) is that narrow spectral lines become smeared into broader line shapes. The 
broadening is a nuisance in some circumstances while in others it provides valued information 
about the distribution of local environments. Interestingly, in some cases quite different 
behavior is predicted to occur. These examples pertain to the quantum states of atoms or 
molecules moving in somewhat heterogeneous environments. Presumably, similar behavior 
may occur for electrons in corresponding situations, e.g. electrons in liquids 

The first example involves atoms or molecules moving in nanoporous media. Nowadays, 
such pores can be nearly perfect and either quasi-one dimensional (ID) and independent 
or interconnected. We focus on the ID case here. The other two examples described in 
this paper involve bundles of carbon nanotubes; these are sometimes called ropes. Such 
bundles consist typically of 100 to 200 tubes, each of which is nearly straight and of diameter 
around 1.5 nm. To see interesting size effects in these problems, the energy levels of motion 
perpendicular to the z axis (the tube or pore axis) should be relatively widely separated. 
That means that the zero-point energy (ZPE), which is of order /i^/(ma^), should be large 
compared to the other energy scales in the problem, such as k^T . Here m is the mass and 
a is the transverse dimension, both of which should be small to reveal the size effects of 
interest. For example, if we discuss H2 in a space with a= 1 nm, the temperature equivalent 
of the ZPE is Ix^ l{mo?kB) ~ \K. 

Section 2 of this paper discusses the problem of particle motion within a nearly cylindrical 
pore. Under appropriate conditions, unusual behavior is found to occur for the adsorbed 
particles' density of states, with implications for the thermodynamic properties of the ad- 
sorbed phase. Section 3 discusses a more complicated case of hydrogen molecules moving 
in the interstitial channels between a triad of nanotubes. Much research has explored these 
states in connection with applications to hydrogen storage and isotope separation, as well 
as intriguing possibilities about ID phases of matter [2-6]. Here, we describe a situation in 
which the heterogeneity associated with an experimental distribution of tube radii yields a 
novel distribution of energy levels, with implications for a predicted bose-einstein condensa- 
tion (BEC) In Section 4, we address briefly a similar problem, the states of molecules 
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strongly bound in the groove region between two nanotubes at the external surface of a 
bundle of nanotubes. Section 5 summarizes our results. 

The evaluation of these examples will, in some cases, employ heuristic models that sim- 
plify the mathematics and clarify the results without sacrificing accuracy of the qualitative 
conclusions. Surprisingly often in science we are lucky and such models turn out to be more 
reliable than they are expected to be. 



2. PARTICLES WITHIN NEARLY CYLINDRICAL PORES 

The pores considered here are assumed to be uniform in the axial (z) direction, so that 
the energy levels of particles moving within them may be written as a sum of a kinetic 
energy of z motion = pi /{2m) of the particles and a transverse energy, Et. 

E^E, + Et (1) 

We first address the transverse problem by finding the energy Ecyi{a) of a particle inside 
a perfectly cylindrical pore (radius a) and then evaluate the effects of eccentricity on the 
spectrum. An atom or molecule confined within a small cylindrical pore experiences a 
potential energy V{r) as a function of radial distance r that can be extremely attractive 
compared to that on a planar surface made of the same material. This expectation is borne 
out in a commonly employed simple model in which the potential is derived by integrating 
Lennard- Jones interactions between the particle and the individual atoms of the host medium 
(r > a). If the parameters characterizing this interaction are e and a, the resulting potential 
energy on the axis of the pore is 

U = V{0) = UoF{R) (2) 

with 

= 1^-15 (3) 

and 

Uo = rnr^e (4) 

Here, n is the number density of the host sohd and R — a/a is the reduced radius of the 
pore. We are especially concerned with low-lying states in this paper. The reduced potential 
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energy function F{R) has its minimum value at R = Rmin = (21/32)^6 ~ 0.932, in which 
case F{Rmin) ~ 0.823. Fig. 1 depicts F{R), which exhibits this minimum, the focus of 
the present discussion. The energy scale for this problem is Uq; if we assume parameters 
corresponding to the H2/graphite interaction {n = 0.11/A'^, a = 2.97A and e = ASK) j^, 
then Uq = 1220i^. This value is about four times as large as the binding energy of H2 on 
a graphite surface and somewhat larger than the low coverage heat of adsorption measured 
on carbon nanotubes. P, 0] 

We proceed to calculate the total energy by solving the Schrodinger equation. In the 
present study, for simplicity, we assume that the potential energy of the particle has a 
constant value V{r) = U for all r < a, but V{r) = 00 for r > a. The ground state 
eigenvalue for this textbook problem is then given by the expression 

Here a = 2.405 is the first zero of the Bessel function of order zero (i.e., the ground 
state wave function). The importance of the second (kinetic energy) term for states near 
R = 1 depends on a dimensionless quantum parameter 77 = h"^ / {nmea^) . For the case of 
the present parameter set, its value is rj = 0.02, implying that the kinetic energy term is 
about 2% of the potential energy, if i? ^ 1. We note that this small value is a consequence 
of the relatively large value of a ~ a. The resulting spread in the wave function is a factor 
about 10 larger than that of H2 in an interstitial channel in a nanotubes bundle, discussed 
in Section 3 below, for which the kinetic energy is a factor ~ 100 greater. 

The total energy function Ecyi{a) has its minimum value for a pore radius near a^m = 
Rmin, i-e. essentially that found in the absence of the kinetic energy term. We call the 
corresponding minimum energy Emin = Ecyi{amin)- We will use an expansion near this 
minimum, 

Ecyl{a) = Emin + ^fJ' {I —f (6) 

Here, /i = (T'^{(PEcyi{a) / da?) with the derivative evaluated at a = amin- Neglecting the tiny 
corrections due to the kinetic energy term, we find /i = 25.6 Uq. 

We now consider the effect of deviation from a perfect cylinder. As an example, we treat 
the case of a quadrupolar distortion, for which the distance R of the pore boundary from 
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the z axis satisfies 



= a[l + Acos(20)] (7) 

(A dipolar deformation would give no energy shift, to this order). Here is the azimuthal 
angle and A is a parameter, assumed small, characterizing the amplitude of the distortion. 
Boundary condition perturbation theory [lO| yields the energy spectrum: 

= E{a, A) = E,yi{a) + (8) 

Here, v = (^a)^/(ma^j„) [1 + a{d In J2{x)/dx)x=a], where J2 is the Bessel function of 
order 2. The numerical value of the quantity in brackets is 1.9. Note that u/fi is of order 
T] <ti I; that is, the energy shift due to the deformation is comparable to the undeformed 
kinetic energy, a small quantity. 

We now consider a particular experimental sample in which the pores' distribution is 
given by a specified function N{a, A). The transverse density of states g{E) for this system 
is derived by integrating over these variables and finding which states have the specified 
energy: 

g{E) = J dadXN{a,X)6[E-E{a,X)] (9) 

We focus here on the lowest-lying states of the problem. These arise from those pores having 
radius near a^m and small values of A. In such a case, we may extract the corresponding 
distribution function as a prefactor: 

g{E) = Niarmn, 0) J da d\ 5[E - E{a, A)] (10) 

The subsequent mathematics is straightforward if one considers very low E, in which case 
one may employ the expansions defined above to perform the integrals analytically: 

giE) = ^iV(a™„, 0)^e{E- E^,„) (11) 



Here 6(x) is the Heaviside unit step function. This equation expresses the fact that the 
transverse density of states is a constant above threshold, at least over the range in which the 
quadratic expansions are valid. This finding may be surprising, at first sight, because one 
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might not have expected a discrete hne spectrum to be broadened into a constant, smooth 
continuum. A phase space argument helps to rationahze the result. The states of a given 
energy E above Emin correspond to the equation fiX^ + = k{E — Emm), where A; is a 
constant and 6=1 —a/amm- This is the equation of an ellipse in the A — 5 parameter space. 
The density of states is proportional to the area between two ellipses of energy differing by 
dE, divided by dE. This analysis leads to the constant density of states specified above. We 
note that this behavior occurs for the same reason that the semiclassical quantum theory 
for the harmonic oscillator yields the correct constant spacing of energy levels, with its two 
degrees of freedom: quantization of the action integral (essentially the area enclosed by the 
ellipse) means that successive states differ by the constant classical frequency. 

We note that the constant (i.e. energy independent) density of states for this problem 
coincides with that of a 2D ideal gas. It is not surprising, therefore, that when this transverse 
spectrum is convoluted with the longitudinal energy contribution E^ the resulting power law 
behavior o{N{E) is that of a 3D gas: N{E) ~ {E-Eminf^- Thus, hydro gen in this confined 
geometry exhibits low energy behavior which it would have in the absence of any confinement 
at all! This conclusion differs from the result for a monodisperse situation, in which case 
the confined particles are characterized by ID physics, i.e. N(E) ~ (i? — Emin)^^^'^ ■ In a 
more complete paper, yet to be written, we will present the implications of this interesting 
result for the thermodynamic properties of the system. There, we will also discuss the case 
of more general deformations of the form cos{p(j)), where p is an integer. It will be seen that 
the quadratic dependence of the energy on the square of the amplitude of the deformation 
is a general property. The effect on the density of states in that general case will differ from 
that described above. 



3. PARTICLES WITHIN INTERSTITIAL CHANNELS 

Bundles of nanotubes and individual tubes are only beginning to be characterized ^| 
and the assignments of chirality (the two integers specifying the kind of tube) and radius 
are usually uncertain. At present, most experimental samples have fairly broad distribu- 
tions of sizes. Nevertheless, most research concerning adsorption has assumed monodisperse 
e„v.o„,„e„ts. for s,mp.,city. An exception to that rule is a reeent paper Q .n wMeh hetero- 
geneity was shown to be essential for understanding adsorption isotherm data [1^ obtained 
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for various gases within bundles. Here, we consider the effects of such heterogeneity on the 
energy spectrum of H2 molecules confined to the interstitial channel (IC) region between the 
tubes. The key function of interest in this analysis is the density of states N{E). As in the 
previous discussion, the assumption of translational invariance in the z direction (accurate 
for nanotubes) means that the total energy of the molecule may be written as a sum of a 
transverse energy {Et) and a longitudinal energy — / (2m), where p is the ID variable 
momentum. 

The interesting problem is the transverse energy, which is very sensitive to the radii (i?i, 
i?2 and R^) of the neighboring tubes. We skip the details of this calculation, which paral- 
lels that of the previous section; a complete description will be published elsewhere. The 
key points to understand are the following. First, there exists a particular combination of 
nanotubes for which the total molecular energy is a global minimum, Emin- This occurs for 
tubes of equal radius {Ri\ — Rmin — 9.95 A, in which case E^nin — —1053K. As in the pre- 
vious section, we first evaluate the transverse density of states g{E) of the molecule. At low 
energy, the relevant states arise from ICs with radii values close to Rmin- To evaluate g{E), 
we determine the iso-energy contours in R space, in which a particular IC is represented by 
a point with coordinates equal to the three radii. Figure 2 depicts the dependence of the 
molecule's energy on a planar surface in R space that passes through the point Rmin{^, 1,1)- 
One observes that these contours are extremely anisotropic; the reason is that the energy 
difference between the global minimum energy and that of the point R = Rmini^, i+x, 1 —x) 
is a quadratic function of x with an extremely small coefficient. This gives rise to the very 
narrow valleys observed in Figure 2. 

In computing g{E), one proceeds to evaluate the number of ICs between energy contours 
at E and E + dE. We find the result that g{E) is proportional to (E — EminY^'^', the power 
1/2 occurs because R space is three-dimensional. (The general behavior in d dimensions is a 
power {d/2 — 1); the previous section exemphfied d — 2). Note that g{E) is mathematically 
coincident with the density of states of a 3D gas. This means that inclusion of motion in the 
z direction gives rise to behavior consistent with an ideal gas in 4D: N{E) is proportional 
to {E — E„iin). Such a spectrum yields an extraordinary consequence: parahydrogen, the 
low temperature (T) form of hydrogen, is therefore predicted to undergo a bose-cinstcin 
condensation (BEC) when confined to this environment, a bundle of nanotubes differing in 
size. The mathematics of this analysis follows that of the textbook derivation of BEC in 



7 



3D.Q| Figure 3 depicts the resulting transition temperature as a function of the molecular 
number density. The typical values are of order milli-Kelvin, easily achievable in a low 
temperature laboratory experiment. 

This predicted transition is remarkable because it is a consequence of the system's het- 
erogeneity. If all of the ICs were identical, in contrast, the description would be that of a 
ID gas, which exhibits no transition. While unusual, this existence of a transition only in 
the inhomogeneous case is not unique. There exist other phase transitions that occur as a 
consequence of nonuniformity, an example of which is the spin glass transition. 

We note several implicit assumptions in the present analysis. One is the neglect of 
intermolecular interactions; this assumption may be valid because of significant screening of 
the intermolecular interactions by the nanotubes An expected singular heat capacity at 
the transition (described elsewhere) relies on the assumption that molecules within an array 
of ICs can redistribute themselves as a function of T. This equilibration may be difficult to 
observe due to the low mobility of the molecules within the ICs. An additional assumption 
is the existence of a continuous distribution of radius values. The radius values of nanotubes 
are determined by the so-called chiral integers (n,m) characterizing the tubes. The spacing 
between consecutive values of radius is less than 0.002 nm, so the discreteness should not 
be a problem. Evidently, the predicted observation^of BEC for hydrogen molecules, not 



seen in any experiment thus far, would be exciting! 
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4. PARTICLES IN THE GROOVE 



The third geometry of interest here is that of the groove region between two nanotubes, 
e.g. at the external surface of a bundle of tubes. We have analyzed this problem by 
diagonalizing the anisotropic dynamical equations of molecular motion in the vicinity of 
the equilibrium position to obtain the two harmonic frequencies {ui and UJ2) characterizing 
this motion, as a function of the two radii. 

Figures 4 and 5 show how these frequencies vary as a function of the radius values. Note 
that the average frequency exhibits a maximum near the point R2 = Ri = 2A. This means 
that the study of the vibrational spectra of these molecules would exhibit behavior charac- 
teristic of a maximum in the vicinity of the frequency corresponding to a double excitation 
(i.e., at Usum = + 002). Without going into detail, this situation is mathematically equiv- 
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alent to that described in Section 2, since that description involves two parameters. Thus, 
the hne shape for an absorption experiment would be constant just below the threshold at 
cOsum, above which no absorption occurs (at least from this combination mode). This is to 
be contrasted with the narrow line present in the absence of heterogeneity. 

5. SUMMARY 

In this paper, we have described a number of inhomogeneous systems that exhibit in- 
teresting physical properties. The common theme discussed here is that the nonuniformity 
leads to anomalous behavior that is qualitatively different from that of the uniform system. 
In some respects, the behavior is much more interesting than that of the uniform system. For 
the physical systems explored here, such heterogeneity is the rule rather than the exception. 
Hence, one should resist the temptation to ignore heterogeneity. 

Joshua Jortner has been an inspiration to one of us (MWC) for many years. We are 
grateful to the National Foundation for support of this research; F.A. acknowledges funding 
from MIUR-COFIN 2001; E.S.H. acknowledges hospitality at the Department of Physics of 
PSU where this work was done. 
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FIG. 2: Contour plot showing the variation Ef^R) — E^in near the minimum Rmin (which is 
chosen as the origin). Prom the contour closest to the minimum (the closed one) to the more 
distant ones, the energy contours correspond to Et(R) — E^in = 0.005,0.01,0.05,0.1,0.5, 1 and 5 
K. 
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FIG. 3: Density of H2 molecules as a function of BEC transition temperature. 
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FIG. 4: Excitation energies (in K) ?iuji, T1U2 of one H2 molecule in the potential well of the groove 
(dashed and dot-dashed lines, respectively) together with the ground state eigenvalue ^(wi +a;2)/2, 
full curve, as functions of radius R (in A) for the symmetric configuration R = R' . 
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FIG. 5: Isoenergy contours of the excitation energies Hui, fiU2 (upper and lower panels, respec- 
tively) as functions of radii i?, R' (in A) . Contours run between ICQ K and 200 K in steps of 4K in 
both graphs. 
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